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Abstract. Models of differentially rotating protoneutron 
stars are calculated, using realistic equations of state of 
dense hot matter. Various conditions within the stellar in- 
terior, corresponding to different stages of protoneutron 
star evolution, are considered. Numerical calculations are 
performed within the approximation of stationary equilib- 
rium, using general relativistic equations of stationary mo- 
tion of differentially rotating, axially symmetric stars and 
using a numerical code based on spectral methods. Fami- 
lies of differentially rotating models of a given baryon mass 
are calculated, using a two-parameter formula describing 
the angular velocity profile within a rotating protoneu- 
tron star. Apart from the usual "mass shedding limit" , we 
introduce an additional "minimal mass limit" for differen- 
tially rotating protoneutron stars resulting from a type II 
supernovae. Maximum angular momentum, which can be 
accommodated by a protoneutron star within these limits 
is calculated, for various thermal conditions in stellar in- 
terior, for a baryon mass of 1.5 M@. In the case of a ther- 
mally homogeneous (isentropic or isothermal) neutrino- 
opaque interior this maximum angular momentum turns 
out to be somewhat higher than that of a cold neutron 
star of the same baryon mass, rotating uniformly at the 
mass shedding angular velocity. However, if the protoneu- 
tron star has a thermal structure characteristic of initial 
state, with a low entropy (unshocked) core, and a high 
entropy (shocked) outer half of baryon mass, the maxi- 
mum angular momentum is significantly lower. This leads 
to a minimum period of uniform rotation of cold neutron 
stars of baryon mass ~ 1.5 Mq, formed directly (i.e. with- 
out a subsequent significant accretion of mass) from pro- 
toneutron stars with shocked envelope, of about 1.7 ms 
and strengthens the hypothesis that millisecond pulsars 
are accretion accelerated neutron stars. 
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I. Introduction 

Newly born neutron stars are expected to be quite dif- 
ferent from ordinary neutron stars, which are a commonly 
accepted model of radio pulsars. Just after their formation 
in a gravitational collapse of a massive stellar core, a future 
neutron star is a hot (internal temperature ~ 10 11 K), lep- 
ton rich (one electron for three nucleons) object - in view 
of these specific features it is called a protoneutron star. 
After some seconds, a protoneutron star transforms into a 
neutron star, loosing the lepton number excess via emis- 
sion of neutrinos trapped in dense, hot stellar interior. The 
short period of neutrino burst is of paramount importance 
for the fate of collapsing star. It is commonly accepted, 
on the basis of numerical simulations of gravitational col- 
lapse of stellar cores, that the prompt-shock mechanism 
fails to make a successful explosion corresponding to a SN 

II. It is currently believed, that neutrino flow from a pro- 
toneutron star can revive the stalled shock. It is thus clear 
that the structure and evolution of a protoneutron star, 
formed as a product of gravitational collapse, are crucial 
for the success (or failure) of the accretion shock revival, 
and therefore, for a successful SN II explosion. 

While the lifetime of a protoneutron star (seconds) is 
negligibly short compared with the lifetime of a neutron 
star at the radio pulsar stage, it is some three orders of 
magnitude longer than the dynamical timescale for these 
objects (milliseconds). In view of this one can study the 
bulk evolution of protoneutron stars in the quasistationary 
approximation, treating convection just as an additional 
transport process in the stellar interior. Static properties 
of protoneutron stars, under various assumptions concern- 
ing their composition and equation of state of hot, dense 
interior, were studied by numerous authors (Burrows & 
Lattimcr 1986, Takatsuka 1995, Bombaci et al. 1995, Bom- 
baci 1996, Prakash et al. 1997, Gondek et al. 1997). 
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The formation of a protoneutron star (PNS) occurs 
on a dynamical timescale. It involves compression, with 
an overshoot of central density, and a hydrodynamical 
bounce. PNS is expected to begin its life pulsating around 
its quasistatic equilibrium. Radial pulsations of PNS were 
recently studied in (Gondek et al. 1997). 

If the presupernova core was rotating - even very 
slowly, the PNS is expected to be rather rapidly rotat- 
ing. This is an inevitable consequence of the conservation 
of the angular momentum of the collapsing core, which 
during a fraction of a second shrinks to about one per- 
cent of its initial radius. The effects of rotation on the 
process of collapse of massive stellar core were studied 
in Monchmeyer & Miiller (1989), Janka & Monchmeyer 
(1989) and Yamada & Sato (1994). 

Rapid rotation can have a strong effect on the structure 
of a PNS. The effect of rotation is stronger than in the case 
of cold neutron stars (NS), because a hot PNS, due to the 
thermal and neutrino trapping effects, is a more extended 
object than a cold configuration of the same baryon (rest) 
mass. Rapid uniform rotation of PNS was studied recently 
in (Hashimoto et al. 1995, Goussard et al. 1997), who con- 
sidered also the consequences of the rotational properties 
of PNS for cold NS, into which a PNS transforms due to 
deleptonization and cooling. 

Uniform rotation of PNS should be considered as a first 
step in the study of their rotational properties. Even if the 
presupernova core was uniformly rotating, collapse would 
generate a significant amount of differential rotation of 
PNS (during collapse, each collapsing annulus conserves 
- to a very good approximation - its angular momentum 
with respect to the rotation axis). Numerical simulations 
of collapse of rotating stellar cores indicate a characteris- 
tic distribution of angular velocity within a newly formed 
PNS : angular velocity of the matter, which is approxi- 
matively constant on coaxial cylinders, decreases with in- 
creasing distance from the rotation axis (Monchmeyer & 
Miiller 1989, Janka & Monchmeyer 1989). 

Evolutionary timescales of a PNS are at most seconds; 
actually, the lifetime of a PNS is shorter than about twenty 
seconds (the time needed for the hot stellar interior to 
deleptonize). The only mechanisms which we could con- 
template as the means of a sufficiently rapid transport of 
angular moment on such a short timescale are related to 
the superstrong magnetic fields, convection and, maybe, 
turbulent viscosity. Our knowledge of the quantitative role 
of these effects, which could rigidify rotation of PNS, is at 
present insufficient to make reliable estimate of rigidifica- 
tion timescale. On the other hand, the viscous timescale, 
related to the transport of neutrinos within the hot, neu- 
trino opaque interior of a PNS, is much longer than the 
PNS lifetime. 

In the present paper we generalize our previous calcu- 
lations of the models of rapidly rotating PNS (Goussard et 
al. 1997) to the case of differential rotation. In particular, 
we study to what extent the presence of differential rota- 



tion modifies the conclusions obtained assuming uniform 
rotation of PNS, referring to the limitations imposed on 
the periods of solitary pulsars formed directly from PNS. 
Existing calculations of rapid differential rotation were re- 
stricted to the case of cold NS (Wilson 1972, Komatsu et 
al. 1989), and used schematic (unrealistic) equations of 
state of T = neutron star matter. Our models of hot, 
differentially rotating PNS are constructed using a realis- 
tic EOS of dense hot matter, taking into account both the 
effects of high temperature, and of trapped neutrinos. 

The model of the dense, hot interior of PNS is de- 
scribed in Sect. 2. We begin with a short description of the 
assumed thermal and lepton structure. Then we present 
the equation of state of the PNS matter, corresponding 
to various stages of the PNS evolution. We discuss also 
characteristic timescales, relevant for PNS, and justify the 
assumption of stationarity (hydrostatic equilibrium) of ro- 
tating PNS. The mathematical formulation of the problem 
of stationary differential rotation of PNS is presented in 
Sect. 3, where we also consider the problem of stability of 
differentially rotating PNS. The numerical method of solu- 
tion of equations of stationary motion is described in Sect. 
4, where we also discuss in some detail the problem of pre- 
cision of our numerical solutions. In Sect. 5 we present our 
models of rapidly, differentially rotating PNS, and discuss 
in particular the mass shedding limits for these models, 
at various stages of PNS evolution. We study also the re- 
lation between differentially rotating PNS and uniformly 
rotating cold NS. Finally, Sect. 6 contains the discussion 
of our results and our conclusions. 

2. Physical conditions in the interior of a pro- 
toneutron star 

Our models of PNS are based on available results of 
(mostly 1-D) simulations of neutron star birth (Burrows & 
Lattimer 1986, Janka & Monchmeyer 1989, Monchmeyer 
& Miiller 1989, Burrows et al. 1995). We consider PNS 
just after it settled into a quasi-stationary state (we ne- 
glect pulsations which could be excited due to hydrody- 
namical overshoot and bounce). As we restrict ourselves 
to the period after the successful revival of the standing 
(accretion) shock, we can neglect the effects of accretion. 
In what follows, we will consider three situations. The first 
one corresponds to a very initial state of a PNS, which is 
then composed of a hot shocked envelope, with entropy per 
baryon (in the units of /cb) s cnv ~ 5 — 10 up to some bary- 
onic density n cnv , and an unshocked core, characterized 
by s coro ~ 1 above some density n coro > n onv . The layer 
with n cnv < n < n coro constitutes a transition region be- 
tween the unshocked core and the shocked envelope. Our 
standard model of a PNS will have baryon (rest) mass of 
Mbar = 1.5 M Q . Then, the shocked envelope will contain 
about 0.7 M Q of the baryon mass. Both core and enve- 
lope are neutrino opaque - a neutrino transparent outer 
layer of the envelope has a negligible mass, and was not 
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included. We will hereafter call such a protoneutron star 
configuration an early-type PNS (EPNS). 

The second situation (established after 0.2-0.5 s, see 
Burrows & Lattimcr 1986, Keil et al. 1996) will corre- 
spond to a thermally homogeneous PNS. At this stage, the 
thermal state of PNS interior (below the neutrinosphere) 
is characterized by s ~ 2, and some characteristic lepton 
fraction Y) . This situation will be referred to as a late- type 
PNS (LPNS). 

Further cooling and deleptonization transforms PNS 
into a NS of the same Mbar- The structure of a sufficiently 
cool NS (internal temperature T < 10 10 K) is determined 
essentially by the T = EOS; due to strong degeneracy 
thermal effects in the EOS of the dense NS interior can 
be completely neglected. This will correspond to a third 
situation considered by us. The three situations, described 
in this paragraph, correspond to quite different EOS of the 
stellar interior, described below. 

2.1. Equations of state and stationary models of PNS 

We will assume a well defined "neutrinosphere" which 
separates a hot, neutrino-opaque interior from colder 
neutrino -transparent envelope. Hot dense matter in the 
hot, neutrino opaque interior is composed of nucleons 
(both free and bound in nuclei) and leptons (electrons, 
positrons, and neutrinos and antineutrinos of all flavors; 
for simplicity we do not include muons). Nuclcon com- 
ponent of hot matter will be described by one of the 
models developed by Lattimer & Swesty (1991), corre- 
sponding to the incompressibility of symmetric nuclear 
matter at saturation density K — 220 MeV. The local 
state of matter in the hot interior is determined by the 
baryon (nucleon) number density n, the net electron frac- 
tion Y e = (n e - — n e +)/n, and the net fraction of trapped 
electron- neutrinos, Y v = (n Vp — n^^/n. An initial stage 
of PNS evolution is characterized by a significant trapped 
lepton number Y\ = Y v + Y e ~ 0.4. Under these condi- 
tions, both electrons and electron-neutrinos are strongly 
degenerate. The deleptonization, due to the v e diffusion 
outward (driven by the gradient of the chemical potential 
of u e ) and subsequent emission from their neutrinosphere, 
implies a strong decrease of Yi on a timescale of seconds 
(Sawyer & Soni 1979, Prakash et al. 1997). After the delep- 
tonization (which for simplicity is assumed to be com- 
plete, with Y v = 0), there is no trapped lepton number, 
so that neutrinos trapped within the hot interior do not 
influence the beta equilibrium of nucleons, electrons and 
positrons. The neutrino diffusion outward (towards the 
neutrinosphere) is then driven by the temperature gradi- 
ent. The neutrinosphere is determined by the condition 
that optical thickness of the layer above it, for electron 
neutrinos with mean energy corresponding to T at the 
neutrinosphere, be unity. 

The effect of high s on the EOS of the shocked envelope 
of a PNS can be seen by comparing the pressure-mass 




11 11.5 12 12.5 13 13.5 



logio P[§ cm " 3 ] 

Fig. 1. Pressure versus matter density in the shocked (s = 5) 
envelope of a PNS (dashed and long-dashed lines). Two curves 
correspond to two limiting cases of trapped lepton number: 
long dashes - maximal trapped lepton fraction Yi — 0.4; short 
dashes - no trapped lepton fraction (Y v = 0). For the sake of 
comparison, we show also the EOS in the standard NS case 
of cold matter (solid line, T = 0), and that corresponding to 
low entropy and no trapped lepton number (dotted line, s = 1, 
Y v = 0). 



density dependence in the low and high entropy cases (Fig. 
1). For s — 5 the pressure increases by more than an order 
of magnitude at p = 10 13 g cm" 3 , as compared to the 
case with s = 1; the relative effect of s on p decreases 
with decreasing density. Above 10 12 g cm" 3 , the effect of 
trapped lepton number on the EOS of shocked envelope 
with s = 5 is much smaller than that of s. 

Actually, our models of EPNS were constructed only 
with a Y v = shocked envelope. This was done mostly for 
the sake of simplicity, and also to counterbalance some 
oversimplifications of our model with constant trapped 
lepton fraction, which would clearly overestimate the role 
of neutrinos (in numerical simulations, Yi decreases out- 
ward, c.f., Keil et al. 1996). However, we think that a 
more detailed treatment of the shocked envelope would 
not change our main conclusions. On one hand, the layer 
with p < 10 12 g cm" 3 , in which the assumption of Y v = 
could lead to a significant underestimating of pressure 
(Fig. 1), contains a very small fraction of the stellar mass. 
On the other hand, the inclusion of Y v > would further 
stiffen the EOS of the shocked envelope, and could only 
make our conclusions concerning the role of the shocked 
envelope even stronger (see Sect. 5.1). 
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2.2. Timescales and stationarity 

The EOS of PNS evolves with time, due to the delep- 
tonization, which changes the composition of the mat- 
ter, and due to temperature change. However, significant 
changes of both T and Yj within the hot, neutrino opaque 
interior occur on a timescale of seconds, if only diffusive 
processes are operating. The evolutionary timescales can 
be shortened, if convective transport is included. Still, 
these evolutionary timescales are orders of magnitude 
longer than the dynamical timescale, governing the read- 
justment of pressure and gravity forces. 

Our treatment of PNS is based on the assumption, 
that convective motions do not influence directly the bulk 
mechanical equilibrium of rotating PNS. For the convec- 
tive motions to be not important dynamically, convective 
matter velocities should be significantly smaller than the 
sound velocity within the neutrino opaque PNS interior. 
Velocity of sound, v s , can be calculated from the EOS of 
the PNS matter using 



7" 



P 



P 



(1) 



where e is the mass-energy density, and the adiabatic in- 
dex 7 is defined by 



7 = 



P 



dp 
drc 



(2) 



Y t ,s 



In the case of cold NS matter the appropriate 7 is that 
calculated at a fixed Y e (Gondek et al. 1997). 

In Fig. 2a we plot the values of v s /c versus matter den- 
sity, within the shocked envelope of a PNS, and compare 
these values with those corresponding to different physical 
conditions in this density interval. The sound velocity in 
the shocked envelope (s = 5) is significantly higher than 
that corresponding to the matter at s = 1. Strongest ef- 
fect of high s is seen for p > I0 12 g cm~ 3 . Trapped lepton 
number increases v g , but its effect is relatively small as 
compared to that of s = 5 for p > 10 12 g cm~ 3 . Sound ve- 
locity in the shocked envelope is significantly higher than 
the velocities of convective motions at the early stage of 
evolution of a protoneutron star obtained in numerical 
calculations (Keil et al. 1996). This conclusion is even 
stronger in the case of the unshocked core (Fig. 2b). The 
effect of finite entropy on the value of v s is there negligi- 
ble, because of the strong degeneracy of the matter, while 
trapped lepton number slightly decreases it. 

A very important timescale is that of viscous dissipa- 
tion of differential rotation. In the absence of convection 
and turbulence, shear viscosity is due to the diffusion of 
neutrinos. We will estimate the neutrino viscosity under 
the conditions prevailing in the dense hot core of PNS, 
using the general scheme of van den Horn & van Weert 
(1981a, b, 1984) for the calculation of transport coefficients 
of neutrinos trapped in dense, hot matter, combined with 
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Fig. 2. a. Velocity of sound (in the units of the velocity of 
light c) versus matter density, in the bloated envelope of PNS 
(dashed line and long-dashed line). The density is bounded 
from above by that corresponding to our assumed value of 
n cnv . The case of Y v = (short dashed line) corresponds to 
no trapped lepton number, while that of Yj = 0.4 (long-dashed 
line) is that of maximum trapped lepton number. For the sake 
of comparison, we show also the case of cold matter (solid line, 
T = 0), and that with low entropy and no trapped lepton 
number (dotted line, s = 1, Y v = 0). 



neutrino mean free paths as quoted in Burrows & Lat- 
timer (1986). We obtain three expressions, corresponding 
to three different physical regimes for dense hot matter 
with trapped neutrinos. The earliest epoch corresponds 
to strongly degenerate electron neutrinos. In the dense 
PNS core, nucleons are degenerate, and neutrino viscos- 
ity in this regime (referred to as Dn, D^), resulting from 
transport of degenerate electron neutrinos, is 



77„(Dn,Di/) 



3.1 • 10 23 P14 x 
x(V70.05)(T/10McV)- 2 



g cm 



-\(3) 



where P14 = p/(10 14 g cm~ 3 ) (our formula is in good 
agreement with the results obtained by Goodwin & 
Pethick (1982), who used the methods of the theory of 
Fermi liquids). After deleptonization, neutrinos become 
non-degenerate, and the neutrino viscosity in degenerate 
nucleon gas becomes somewhat lower, 



?7, y (Dn,NDi/) ~ 9.7- 10 22 x 

x(T/10 MeV)/pif g cm" 1 s" 



(4) 



Finally, at subnuclear densities we can have non-degene- 
rate neutrinos trapped in a dense hot nucleon gas, which 



J-O. Goussard et al.: Rapid differential rotation of protoneutron stars and constraints on radio pulsars periods 



5 



0.8 



.0.6 



0.4 




I 

10 



15 

p[10 14 g cm" 3 ] 



20 



25 



Fig. 2. b. Velocity of sound (in the units of velocity of light 
c) versus matter density, within the protoneutron star core 
(n > n COI . c ). Notation as in Fig. 2a. 



corresponds to 
77„(NDn,NDi/) 



2.3 • 10 23 x 
x(T/10 MeV) 2 /pi3 g cm" 



(5) 



In all cases the viscous timescales within the stellar inte- 
rior, r„ = I 2 p/riv, with / ~ 10 km, are much longer than 
the lifetime of a PNS. Both magnetic field effects and lo- 
cal turbulence could shorten viscous timescales, but they 
are expected to be still much longer than the dynamical 
timescale. 

The above discussion shows, that evolutionary time- 
scales within a PNS, resulting from the transport pro- 
cesses in the stellar interior, are significantly longer than 
the dynamical timescale for PNS. In view of this, we are 
able to decouple the PNS evolution from its dynamics. We 
can then treat PNS rotation in a stationary approxima- 
tion, with a well defined EOS of the PNS matter. 



3. Formulation of the problem 

The problem of the calculation of the stationary state of 
differentially rotating cold neutron stars, within the frame- 
work of general relativity, was considered previously by 
some authors (Komatsu et al. 1989, Wilson 1972, 1973), 
who restricted themselves, however, to polytropic (and 
hence unrealistic) EOS of dense matter. Here we will ex- 
tend the methods used at T = to the case of hot PNS 
with realistic EOS. We will use the notation and formalism 



developed in Bonazzola et al. (1993), hereafter referred to 
as BGSM. 

3.1. Equation of stationary motion 

As we already mentioned in paper I, one of the problems 
introduced by finite temperature and differential rotation, 
if one does not make any further assumption on the equi- 
librium configuration, is that the equation of stationary 
motion does not have a first integral (Bardeen 1972). This 
equation reads (eq. (5.8) of BGSM) : 



dip 
e+p 



cUn( 



r ' 



-Fdifl 



(G) 



where p is the pressure, e the mass-energy density, N the 
lapse function appearing in the space-time metric, T is the 
Lorentz factor due to rotation, £1 is the angular velocity 
and F = u$u l where u is the fluid 4- velocity. 

This equation is the general relativistic equivalent of a 
well-known Newtonian formula (see e.g. Tassoul 1978). It 
is straightforward to show that a sufficient condition for 
(p|) to be integrable is 



T = T(n) . 

Eq.(^|) can then be rewritten 

d i [G + ln(^)] = -Fd i Q, 

where the heat function G is defined by 
dp 



G = 



i(n,T(n)) + p(n,T(n)) 



(7) 



(8) 



(9) 



which is not to be confused in non-isentropic cases with 
the logarithm of the dimensionless relativistic enthalpy 
per baryon H = ln(/) = ln[(e + p) / (nmoc 2 )} (see paper I). 
Then, as in BGSM (§5.1), we apply the Schwarz theorem 
and find that cither 



Q, = const. 



F = F(Q). 



(10) 



(11) 



This corresponds to the well-known result that, assuming 
barotropic equilibrium (T = T(n)), the rotation of the 
star is either uniform or cylindrical. Being interested in 
differential rotation, we restrict ourselves to the case (pi]). 

The rotation law is obtained by solving Eq. (5.17) of 
BGSM (we correct here a misprint of this paper) : 



F(fi)- 



A 4 B 2 r 2 sin 2 6(fl - N<t>) 
N 2 - A 4 B 2 r 2 sin 2 0(fi - N<t>)'< 



= 0, 



(12) 



where A, B are metric potentials and is the shift vector 
(see BGSM for a precise definition). 
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One then obtains the first integral of motion in the 
form : 



N f 

G + Ln(— ) + / F(Q)dn = const. 



(13) 



This enables us to calculate G in the whole star and 
thus the other thermodynamic quantities, using the G- 
tabulated EOS of the form e = e(G), p = p(G), s = s(G) 
(see §4.2). 

Finally, in the case of general relativistic thermal equi- 
librium, one can show, in the context of the Eckart (1940) 
or Carter (1983) theory of heat diffusion, that 



T* = T-ef Fdn = const. 

r 



(14) 



where T* is the redshifted temperature which would be 
measured by a distant observer, corresponds to thermal 
equilibrium (no heat flux). Such a T-profile leads to a 
simple form for the first integral of motion : 



const. 



(15) 



where pL is the baryon chemical potential. The T-profile is 
then determined from 



/j,(n,T(n)) 
T(n) 



const. 



(16) 



which allows us to use the previous formulation (|l^) of 
the first integral of motion. 

3.2. Differential rotation law 

There is obviously an infinite number of degrees of free- 
dom in the choice of the F-l&w. However, we wish to have 
a law which: (1) approximately reproduces the f2-profile 
obtained in 2D-simulations, i.e. essentially decreases to- 
ward the edge of the star; (2) is a "simple" law, in order 
to limit the number of free parameters. We thus chose the 
rotation law used by Komatsu et al. (1989) : 



F(n) = Rl(n c - n) 



(17) 



where O c is the central angular velocity and i? 2 , is a free 
parameter, with the dimension of the square of a length 
(denoted A in Komatsu et al. 1989). This simple rotation 
law allows for a rapid solution of (|lj) since it becomes 
then a cubic equation. 



In the Newtonian limit, the rotation law (17) yields 
the angular velocity 



Rl 



r 2 sin 2 i 



(18) 



We see that Rq is thus a characteristic scale of spatial 
variation of Q : when Rq is small the rotation profile is 



very steep, leading to a constant specific angular momen- 
tum distribution, and when it is large the rotation is al- 
most uniform. Numerical experiments have shown that the 
general relativistic profile has the same qualitative prop- 
erties as the Newtonian one. Note that this Newtonian 
profile is also the equivalent, for cylindrical rotation, to 
the one used by Monchmeyer & Muller (1989) and Janka 
& Monchmeyer (1989) at the beginning of their calcula- 
tions with spherical rotation. It is also the one used by 
Muller and Eriguchi (1985) in their calculations of "fiz- 
zling" white dwarfs. 

3.3. Stability of differentially rotating PNS 

Differentially rotating self-gravitating bodies are subject 
to various dynamical instabilities, which may deeply affect 
their structure (see e.g. Zahn 1993). In what follows, we 
study to what extent our models of hot PNS are stable re- 
garding the two main axisymmetric instabilities due to the 
differential character of rotation, appearing in barotropic 
configurations : the Rayleigh and the shear instability. We 
then briefly discuss the case of secular instabilities linked 
to the spontaneous symmetry breaking phenomenon. 

3.3.1. Rayleigh instability 

It is straightforward to show that the newtonian ll-profile 
Lq) satisfies the classical Rayleigh stability criterion, i.e. 



_ d(r 2 sin 2 0fi) 
&(rsm0) d(rsin0) 



>0, 



(19) 



where j is the specific angular momentum of baryonic mat- 
ter. This criterion yields also a sufficient condition for sta- 
bility in the case of cylindrical rotation laws, i.e. 



^<0 



(20) 



and one can show that the general relativistic Sl-profile 
obtained with the law ([l7]) fulfils this last condition in the 
isentropic case (Komatsu et al., 1989, §3.1). For the full 
general relativistic formulation in non-iscntropic cases, we 
checked that all our models verified (|lj]) . As this stability 
analysis is a local one in which the gravity does not enter, 
we infer that ( fl9| ) is a good criterion even in the general 
relativistic case. 

3.3.2. Shear instability 

The case of the shear instability is much more compli- 
cated, and we will restrict ourselves to crude estimates. 
Our configurations, in which the Reynolds number (which 
compares turbulent to molecular viscosity) is of the order 
of Re = R 2 pfl/n ~ 10 3 (see Section 2.2), are very likely to 
be unstable (see Zahn 1993). In principle, this instability 
could be blocked by a stable entropy or specific angular 
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momentum stratification inside the star. However, even if 
the entropy and specific angular momentum gradients of 
our models were all stable (apart from a thin outer layer 
of the star), these gradients are too weak to stabilize the 
flow. Indeed, the Richardson number Ri, which compares 
the efficiency of the stabilizing s and j-gradients to desta- 
bilizing effect of shear, can be written in the Newtonian 
case as 



Ri 



9- — )\r 

dU 

"d7 



ldln(p) dln(p) 
7 dr dr 



1^! 

r 3 dr 



(21) 



where g is the gravitational acceleration, 7 the adiabatic 
index, p the mass density and dU/dr the vorticity. In view 
of the fact, that our models are almost adiabatic, i.e. that 



ldln(p) dln(p) 
7 dr dr 



< 0.5 



(22) 



in the whole star, and that as i?g decreases the j- 
distribution tends to a j = const, one (giving thus a 
marginally Rayleigh-stable stratification) , it was not a sur- 
prise to find that almost all our models verified Ri < 
Riait — \ (at least in some part of the star), which is 
a sufficient condition for instability. 

However, the global effect of the shear instability is to 
suppress the phenomenon that causes it, namely, differ- 
ential rotation. As the maximum angular momentum of 
PNS decreases with decreasing differential rotation (see 
§5), the introduction of this instability can only decrease 
the maximum angular velocity of such stars. 

3.4- Non- axisymmetric instabilities 

A rapidly rotating PNS can be affected by secular instabil- 
ities, driven by the gravitational radiation reaction (GRR 
instability) or by viscous dissipation. However, detailed 
calculations performed for hot uniformly rotating NS sug- 
gest that, for temperatures exceeding 10 10 K, the GRR 
instability can only slightly decrease the maximum rota- 
tion frequency (Cutler et al. 1990, Ipser & Lindblom 1991, 
Lindblom 1995, Yoshida & Eriguchi 1995, Zdunik 1996). If 
the instability is driven by viscous dissipation, Bonazzola 
et al. (1995) have shown for cold uniformly rotating NS 
that only very stiff equations of state are likely to allow 
for spontaneous symmetry breaking. 

As our models are differentially rotating and hot, no 
previous study applies to our EPNS and LPNS. To es- 
timate to what extent our models could be plagued by 
spontaneous symmetry breaking, we used the newtonian 
criterion of stability T/|VF < 0.14, where T is the ki- 
netic energy and W the gravitational energy of the star. 
As our PNS models are slowly rotating compared to cold 
NS, we were not surprised to find that T/|W| <C 0.14 for 



the EPNS models. For LPNS, we found some models for 
which T/|W| ~ 0.12 — 0.14, which, considering the approx- 
imate criterion we used, cannot be viewed as conclusive of 
symmetry breaking. Furthermore, inclusion of symmetry 
breaking effects can only decrease the maximum angular 
momentum of PNS and thus, their maximal rotational ve- 
locity, which would only strengthen our conclusions. 

4. Numerical method 

4-1. Axisymmetric code 

We used the same code based on spectral methods as in 
paper I, modified to take into account differential rotation. 
We refer the reader to our previous paper for a brief de- 
scription of the code and tests of precision, which we made 
in the case of uniform rotation. Here we briefly outline the 
different or new features of the code. 

The global relative error, evaluated by the means of 
the virial identity test (see BGSM), is ~ 10~ 5 . This is by 
two orders of magnitude better than the precision reached 
in paper I. Such an increase of precision was possible due 
to new interpolation procedure used for the EOS (see next 
section). 

We used two grids in r for the star and one in the 
1/r- variable for the exterior, of N r Chebychev coefficients, 
with N r — 65 in the interior and N r = 33 in the exterior, 
and one grid in 9 with Ng = 33. We checked that taking 
a greater N r or Ng does not change the results by more 
than a few 10~ 6 at most, which stays within the global 
precision reached. Let us stress that such a low number 
of grid points is sufficient to reach high accuracy within a 
spectral method (this would not be possible in the case of 
a finite difference scheme). 

The temperature drop at the "neutrinosphere" is not 
always located on the border of the internal grid in r, 
which could influence the precision (remember that the 
spectral methods are very sensitive to discontinuities). We 
checked that, in fact, even when the "neutrinosphere" is 
far from the border of the grid, the global physical proper- 
ties of the star did not change much for the late-type mod- 
els of PNS (we found also relative variations of a few 10~ 6 
at most). For early-type models of PNS, the results were 
more sensitive to the location of the "neutrinosphere" , and 
we were obliged to position the grids quite precisely to 
reach the former precision. 

4-2. Interpolation procedure for the EOS 

The precision of the axisymmetric code is very sensi- 
tive to the thermodynamic consistency of the interpolated 
EOS. In particular, the following thcrmodynamical rela- 
tion must be fulfilled to ensure a good precision of the 
code, 



dp 

dG =£ + p - 



(23) 
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If the above relation is only approximately fulfilled by the 
EOS (as is the case when various thermodynamical quan- 
tities are independently interpolated) the equation of sta- 
tionary motion cannot be strictly satisfied, neither. This 
influences greatly the value of the virial parameter and 
thus, the precision of the calculations. In order to have a 
consistent EOS (in the sense that §2^) is verified), we used 
an original procedure, adapted from a method derived by 
Swesty (1996). 

Once the temperature profile is chosen, the EOS can 
be tabulated. The heat function G is calculated by inter- 
polating the quantity l/(e + p) in the integrand in Eq. (4) 
as a function of p using cubic splines and integrating this 
interpolating polynomial. The quantity G can then be cho- 
sen as the new interpolation variable : as this is a quantity 
given directly by the first integral of motion (see §3) this 
ensures a rapid calculation of the thermodynamical quan- 
tities. Futhermore, following Swesty, we use the fact that 
dp/dG is known and interpolate p and dp/dG by means 
of a Hermite interpolation scheme of degree 3. The energy 
density is obtained at each point by derivating the inter- 
polant and using (^3|) . Other thermodynamical quantities 
can be interpolated independently since they do not enter 
this relation. In doing so, we ensure that ( p3| ) is valid for 
the whole interpolation polynomial and not only at the 
interpolation points. 

With such an interpolation procedure, we were able to 
lower the virial error parameter by two orders of magni- 
tude. 

5. Numerical results 

5.1. Static models of early-type PNS and "minimal mass" 
limit 

Early-type PNS, with a high-entropy, bloated envelope, 
are short living structures. Their lifetime might be as short 
as 100-200 ms (if one includes convective motions, see Keil 
et al. 1996), but this is still a long time on a dynamical 
timescale. As early-type models of PNS have not been in- 
vestigated in paper I, we study here the static properties 
of these objets. The gravitational and baryonic mass of 
static models versus central baryonic density are displayed 
in Fig. ||. Let us notice, that we have restricted ourselves 
to the range of masses, consistent with assumed thermal 
structure of PNS. As one can see, there is a minimal mass 
for such equilibria, as for cold ones, but it is of the order 
of M™(EPNS) - 1.4 M Q , much higher than the stan- 
dard minimal mass of cold NS (M™(NS) - 0.1 M Q ). 
This high minimal mass of EPNS is due to the fact that 
thermal pressure at the base of the neutrinosphere tends 
to blow-off the envelope of the star, and thus additional 
mass compared to the cold case is needed to preserve equi- 
librium. This effect was already present in the models of 
paper I, but we were at that time interested in the max- 
imal mass of PNS and did not investigate this problem. 
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Fig. 3. Baryon and gravitational mass of static models of 
early- type PNS versus central density, n ccntr . The solid curve 
is for baryonic mass and the dashed curve for gravitational 
mass. The upper long-dash/short-dash curve is the baryon 
mass, Afbar , calculated for differentially rotating configurations 
at Rq = 1 km 2 and RQ C — 0.284 c, where R is the star radius. 
This illustrates the effect of rotation on the minimal bary- 
onic mass. The parameters of the EOS chosen are s cnv = 5, 
Score = 1 ka, n cnv = 0.02 fm~ 3 and n coro = 0.26 fm~ 3 , with 
Yi — 0.4 in the core and Y v = in the shocked envelope. 

For late- type PNS, the minimal mass is increased at most 
up to M™ n (LPNS) - 0.9 M Q (Gondek et al., 1997). 

The value of M^; n (EPNS) is strongly linked to the 
properties of our EOS, through the assumed values of the 
EOS parameters. Table 1 displays the minimal masses 
reachable for various combinations of the parameters of 
the EOS. 

For some sets of parameters, minimal baryonic mass 
is larger than 1.5 Mq. In particular, we found that for 
s cnv > 6, the minimal baryonic mass is greater than 
1.5 M© and therefore no static equilibrium can be found 
for our "canonical" baryon mass. Let us notice, that the 
value of minimal mass is not sensitive to the presence of 
trapped neutrinos in the unshocked core : had we taken 
a neutrino free core, the minimal mass would change only 
by a few percent as compared to that calculated assuming 
Y { = 0.4. 

Let us stress, that EPNS models with n cen tr lower than 
that corresponding to the minimum mass configuration 
(i.e., those to the left of the minima of the mass - central 
baryon density curves in Fig. 3) are rather unrealistic, and 
moreover likely to be unstable. They consist mostly of the 
shocked, high entropy matter. Clearly, such objects could 
not be formed in a successful SN II explosion: too much 
of the shock energy would be dissipated for heating and 
dissociating the collapsing matter. 

The very existence of a high minimal baryonic mass 
for static configurations of early-type PNS leads to severe 
constraints on their rotational velocity or, more precisely, 
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Table 1. Minimal baryonic masses of EPNS for various values of the parameters of the EOS. 



Scnv 5 4 

Score 1 1 

n cnv [fm" 3 ] 0.02 0.02 

n C ore[fm- 3 ] 0.26 0.26 
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0.02 


0.02 


0.02 


0.01 


0.26 


0.3 


0.1 


0.26 



M bar , min [M ] 1.41 1.23 1.56 1.50 1.19 1.22 



on the maximal angular momentum they can bear. As one 
increases the central angular velocity of such stars, keeping 
their baryonic mass constant, the centrifugal force adds 
to the thermal pressure at the base of the envelope. The 
envelope is then more easily lifted than in the static case, 
leading to an increase of the minimal mass of stationary 
configurations. At the "minimal mass" limit, the minimal 
mass coincides with the canonical one and no equilibrium 
can be found with a higher angular momentum. So, apart 
from the usual keplerian - or "mass shedding" - limit (for 
which the centrifugal force exactly balances the gravity 
at the equator), this "minimal mass" limit also constrains 
the domains of variations of il c and Rq. 

5.2. Variations of J versus O c for early-type PNS 

Given the two limits constraining Jl c and i?Q, the prob- 
lem remains of finding the maximal angular momentum of 
EPNS. Variations of the angular momentum versus cen- 
tral rotational velocity were previously studied for cold 
uniformly rotating NS with various realistic EOS (Cook et 
al. 1992, 1994, Salgado et al. 1994). These studies showed 
that there exist two distinct families of uniformly rotat- 
ing cold NS : normal and supra-massive stars. Normal 
stars have a baryonic mass lower than the maximal static 
baryonic mass, and their angular momentum is an increas- 
ing function of the central angular velocity. Supra- massive 
stars have baryonic mass greater than the maximal static 
baryonic mass, which they would thus not support with- 
out rotation. Their angular momentum is not a monotonic 
function of the central angular velocity, and the stable ke- 
plerian configuration is not that of maximal angular mo- 
mentum. This property leads so to the well-known phe- 
nomenon of "spin-up by angular momentum loss" . 

As the canonical mass we chose (Mb ar = 1-5 M©) is 
less than the maximal static baryonic mass for all the 
EOS we used, we expected that, for a given Jl-contrast (a 
given Rq), J would increase with increasing central angu- 
lar velocity for our models of PNS. Numerical experiments 
showed that this is indeed the case for cold models but not 
for the hot ones. This can be understood as follows. For 
cold uniformly rotating NS, the phenomenon of spin-up by 
angular momentum loss is due to the high oblateness of the 
supra-massive stars due to their high rotational velocities 
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Fig. 4. a. Angular momentum versus central angular velocity 
of EPNS of baryon mass 1.5 Mq, for rigid rotation and for 
differential rotation with Rl = 1.2, 1.1, 1.0, 0.9, 0.8 and 0.7 
km 2 (from the left to the right). The parameters of the EOS 
are as in Fig. fa. Triangles denote the keplerian configurations. 



(see Cook et al., 1992, 1994) : since 51 = J/I, the varia- 
tions of SI depend on the variations of both J and I, for 
a sufficiently extended star, the redistribution of / when 
one increases 51 can lead to a decreasing J. This effect is 
thus strongly linked to the radial distribution of matter 
inside the star, and can be significantly enhanced by the 
presence of differential rotation (Shapiro et al., 1990), but 
nevertheless appears only in supra-massives stars. This is 
not the case for PNS. Due to their extended structure in- 
duced by the supplementary thermal pressure, very little 
rotation is needed to produce this phenomenon. Thus it 
can appear even in normal PNS. 

Figures |^a and ^b depict typical behaviours of J as a 
function of Sl c : starting from the static model, we increase 
the quantity i?Sl c (which parametrizes our models) up to 
the limiting configurations, keeping the baryonic mass and 
i?§ constant. As one can see, the "spin-up" phenomenon 
is not present for uniformly rotating shocked PNS. Dif- 
ferential rotation has the effect of increasing the maximal 
central angular velocity and the maximal angular momen- 
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Fig. 4. b. Angular momentum versus central angular veloc- 
ity of early- type PNS, of baryon mass 1.5 M©, for differential 
rotation with Rq = 0.6, 0.5, 0.4, 0.3 and 0.2 km 2 (from the 
left to the right). The parameters of the EOS are as in Fig. 
^| Triangles denote the keplerian configurations and circles are 
for configurations with M™™ = 1.5 Mq ("minimal mass" limit 
configurations). The upper-left curves are the J ms < J < J mm 
parts of the i? 2 = 0.6, 0.5, 0.4 km 2 curves (see text). 



turn J max , "pushing" the curves to the right and the top 
in Fig. [| It also increases the angular momentum J ms 
at which mass shedding appears (configurations denoted 
by triangles on Fig. 0a) : J ms is an increasing function of 
l/Rf, 

The two physical limits previously mentioned restrain 
the available maximal angular momentum. For small Q- 
contrast (large Rq), the mass shedding limit is reached 
first, and one cannot store enough angular momentum in 
the star to reach the "minimal mass" limit at J = J mm 
: this is the case for all the configurations of Fig. ^a, for 
which the angular momentum J verifies < J < J max , 
with J ms < J max < J mm . Let us stress that in the case 
of cold uniformly rotating normal NS, one has always 
jmax _ £ or a fi xec [ Mb ar , which is never the case for 
the differentially rotating sequences displayed in Fig. 4a, 
b. 

As we decrease Rq, J max increases and the J max con- 
figuration will finally reach the minimal mass limit for 
J max = J mm ~lM kmc for the EOS we chose. With 
a further increase of the fi-contrast (decrease of Rq), the 
J — f2 c plot at fixed Mb ar divides into two disconnected 
segments (branches). The first branch, which starts at 
r2 c = 0, ends at the minimum mass limit with J mm . A 
second segment (upper left curves in Fig. 0Jd) is bounded 
from above by another minimum mass limit, </ m m (which 
is actually an overall maximum of J for a specific value 
of -Rq), and from below by the mass shedding limit with 
J ms . Notice, that these segments evolved from the upper 
left parts of the J-f2 c curves in Fig. 0a, with the very 



upper part removed by the minimum mass constraint. 
The upper left segments correspond to rotating config- 
urations with J ms < J < J,' nm = J max . Such a situation 
is characteristic of a rather narrow range of the Rq pa- 
rameter (Rq — 0.6, 0.5 and 0.4 km 2 in Fig. ||d). For a 
still larger fi-contrast (still smaller Rq) the left segment 
disappears and we are left only with the configurations 
verifying < J < J mm = J max (curves with i? 2 , = 0.3, 0.2 
km 2 in Fig. 0b). 

We show in Fig. |5|a and ||b the constant energy density 
lines in the cross sections of two typical limiting configu- 
rations of EPNS. As one can see, the keplerian configu- 
ration is a very extended one, with a "disk-like" shape 
due to the combined effect of differential rotation and of 
a high entropy envelope. The minimal mass limit configu- 
ration is more compact, with a "toroidal" shape. In both 
cases, small irregularities in the outermost (low density) 
constant energy density lines are due to the Gibbs phe- 
nomenon in the 8 direction, a problem inherent in our 
numerical scheme but small enough not to affect badly 
the virial parameter (see sect. 4.1). 
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Fig. 5. a. Iso-energy density lines of the keplerian configura- 
tion of Mbar = 1-5 Mq, with Rq = 0.5 km 2 (corresponding 
to the triangle of the Rq = 0.5 curve in figure ^b). x and y 
are related to the MSQI coordinates r and 6 by x = rcosO, 
y = rsinfl (see BGSM for details). 



5.3. Thermal structure and maximal angular momentum 
of PNS 

Uniformly rotating models of late- type PNS, considered 
in Paper I, assumed a relatively homogeneous distribu- 
tion of entropy within the neutrino opaque interior. In 
order to answer the question, to what extent the conclu- 
sions of Paper I would be modified by allowing for dif- 
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EOS s corc = 2 s corc = 1 T* — 15 MeV T* = 25 MeV 

Yi = 0.4 Y =0.4 Y u = Y v = 



Jig crit [km 2 ] 1.044 1.091 
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Fig. 5. b. Iso-energy density lines of the minimal mass limit 
configuration of Mb ar = 1-5 Mq, with Rq — 0.5 km 2 (corre- 
sponding to the upper right circle of the Rq — 0.5 curve in 
figure |b). 



ferential rotation, we considered models, which were isen- 
tropic or isothermal within the neutrino opaque core. Let 
us mention, that such a homogeneous thermal state of a 
PNS corresponds to a relatively late switching-off of accre- 
tion due to the lift-off of the shock. Accretion stage had 
to be sufficiently long to allow for settling-down of the 
PNS envelope, and for the thermal uniformization of the 
neutrino-opaque interior. In what follows, we will restrict 
ourselves to the case of Afb ar = 1-5 M Q . As in the case of 
two-component EPNS, described in the preceding section, 
we calculated J(f2 c ) sequences for fixed parameter Rq. We 
found, that there exist a specific value, i?o,crit, such that 
for Rq < i?o,crit we get the maximum angular momentum 
of PNS which is larger than the maximum angular momen- 
tum of cold, uniformly rotating NS of the same Mbar- So, 
for Rq < -Ro.crit, differential rotation allows PNS to have 
enough angular momentum to yield cold NS, rotating uni- 
formly at maximal angular velocity. The values of i?o,crit 
for the different EOS used in paper I are displayed in table 
2. It should be stressed, that this qualitative conclusion is 
generally valid for thermally homogeneous LPNS, for all 
considered EOS of stellar interior (isentropic with s < 2, 



isothermal, with no trapped lepton number Y v = 0, and 
with maximum trapped lepton fraction Yi — 0.4), and 
that the actual value of i?o, C rit is roughly independent of 
the assumed EOS. 

For Rq < -Ro.crit ! we can consider configurations having 
an angular momentum greater than the maximal angular 
momentum of cold NS. Such PNS will necessarily eject a 
part of J and Mbar during their evolution toward cold, 
uniformly rotating NS. This will probably lead to a max- 
imally rotating (i.e., rotating at its mass shedding limit) 
NS surrounded by an ejected matter. For Rq > i?o,crit, 
evolution of PNS into a cold NS leads to the final rotation 
frequency which is lower than the maximal one for such a 
cold NS, but the difference decreases with decreasing Rq. 

Let us consider now the case of early- type PNS. Be- 
ing aware of the particular behaviour of J depicted in 
the previous paragraph, we scanned the whole domain of 
0.1 < Rq < 1 km 2 for the maximal angular momentum 
of LPNS. In this domain, the maximal angular momen- 
tum configurations were never keplerian. We found that 
1.0 < J max < 1.16 M Q km c, which is roughly indepen- 
dent of Rq. As the crucial factor determining the "mini- 
mal mass limit" is the effective gravity at the base of the 
shocked envelope for large Jl-contrast, it is not surprising 
to find that our results are roughly independent of the pa- 
rameter Rq (thus, in some way, of the assumed rotation 
law), and depends only on the angular momentum of the 
star. These limits on the total angular momentum of the 
early-type PNS leads to limits on the period of cold NS of 
the same baryon mass ofl.6<P<1.8 ms. 

These results are relevant for NS which originate in a 
type II supernova and evolve without significant accretion 
from a hot, differentially rotating LPNS to a cold, uni- 
formly rotating NS. If the baryonic mass of such a PNS is 
~ 1.5 Mq, then the maximum frequency of the resulting 
NS cannot exceed ~ 600 Hz. This evolutionary argument 
thus imposes more stringent limits on the rotation rate of 
NS than standard calculations (Cook et al. 1994, Salgado 
et al. 1994). 

6. Discussion and conclusions 

In the present paper we studied models of rapidly differ- 
entially rotating protoneutron stars. Using realistic equa- 
tions of state of the protoneutron star interior, we cal- 
culated, using general relativistic equations of stationary 
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motion, families of configurations corresponding to the 
main stages of evolution of protoneutron stars. The ear- 
liest stage corresponds to two-component objects, with a 
low entropy unshocked core and a bloated, high-entropy 
envelope. The next stage is that of a lepton rich protoneu- 
tron star, with a relatively uniform distribution of entropy 
per baryon within the hot neutrino-opaque core, contain- 
ing degenerate neutrinos. Finally, the latest stage is that 
a deleptonized hot neutrino-opaque core, which contains 
non-degenerate neutrinos. We considered mostly the case 
of a protoneutron star of the "canonical baryon mass" 
1.5 Mq, which roughly corresponds to measured gravita- 
tional masses of binary pulsars. 

Configurations of stationary motion were calculated 
assuming axial symmetry, and neglecting effects of merid- 
ional circulation and convection. Such an approximation 
was valid, because the velocities of meridional circulation, 
as well as those connected with convective motion, were 
much smaller than the velocity of sound in the protoneu- 
tron star interior. On the other hand, they were also much 
smaller than the linear equatorial velocities considered. 
While convective motions could be assumed to have no 
direct dynamical effect, they influenced the PNS structure 
through changes in composition and temperature, through 
transport of heat and of lepton number, leading to the 
evolution of the equation of state of stellar interior. Our 
assumption that the local temperature was a function of 
local density was sufficient for the equations of stationary 
motion to be integrable; stationary state of differential ro- 
tation was then that of a barotropic equilibrium. 

The approximation of stationarity was crucial for the 
simplicity of calculations. The shear viscosity, resulting 
from neutrino diffusion, implies angular momentum trans- 
port and heating of the stellar interior. However, the 
estimates of the characteristic timescales, under various 
regimes prevailing within the hot neutrino opaque core, 
imply that these effects can be neglected on the timescale 
of the PNS lifetime. Clearly, the timescale of angular mo- 
mentum and energy transport can actually be significantly 
shortened by convective and turbulent motions within the 
stellar interior. Also, the presence of magnetic field, and 
subsequent amplifying of it due to the differential charac- 
ter of rotation (winding-up of magnetic field lines), could 
contribute to the angular momentum transport and to 
rigidifying of the rotational motion. However, even in this 
case characteristic timescales referring to the global struc- 
ture of rotating star may be expected to be of the order 
of seconds. Clearly, all these effects which are expected to 
be operating in a differentially rotating protoneutron star, 
deserve a separate study. 

Our models describe rotating PNS after revival and 
lift-off of the shock, i.e., after accretion onto PNS became 
insignificant for mechanical equilibrium. Basically, we con- 
sidered two limiting scenarios of the 'liberation' of PNS 
from accretion. The first case corresponds to a relatively 
early switching-off of accretion, so that rotating PNS has 



still a pronounced two-component thermal structure, with 
an unshocked low entropy core and a shocked high entropy 
envelope. The second case is that of a relatively late lifting- 
off of the shock, when the accretion stage was sufficiently 
long to allow for mixing of the low and high entropy mat- 
ter, leading to a roughly isentropic PNS interior. Let us 
notice, however, that the second scenario is, in view of 
very high accretion rate, quite risky for the PNS survival. 
Huge accretion rates (of the order of M /s) might result 
in a high probability of exceeding the maximum allowable 
mass for neutron stars, and therefore could easily lead to 
formation of a black hole. 

For a "canonical baryon mass" of 1.5 M , and within 
the considered family of differentially rotating configura- 
tions, we find a set of configurations with angular momen- 
tum exceeding the maximum one that can be accommo- 
dated by rigidly rotating cold neutron stars of the same 
baryon mass. We expect, that the evolution of such differ- 
entially rotating protoneutron stars into cold, uniformly 
rotating neutron stars will be accompanied by ejection of 
matter from their equators. Such a mass shedding config- 
uration may be expected to be susceptible to some addi- 
tional instabilities, which - if breaking axial symmetry - 
could contribute to a further angular momentum loss via 
emission of gravitational radiation. 

Assuming that the initial differentially rotating con- 
figuration had a characteristic two-component structure, 
with a high entropy envelope, containing some half of the 
baryon mass of the protoneutron star, we get (for a stan- 
dard baryon mass of 1.5 M Q ) a minimum period of result- 
ing cold neutron star P^J-Jf — 1.7 ms. Shorter rotation pe- 
riods of cold neutron star of the baryon mass of ~ 1.5 Mq 
(such as the periods measured for the most rapid mil- 
lisecond pulsars like PSR 0034-0534, PSR 1937+21, PSR 
1957+20, see e.g. Bailes & Lorimer 1995) could then be 
reached only by a spin-up via accretion of a significant 
fraction of solar mass from an accretion disk. 
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